function [outputArg1] = mkfigsix(xdates,spotprices,fdates,futures,modelfut,modelfutci,othermodelfut,legendtext); 
%mkfigsixbw Make the plot for LMV RESTAT.
%   Given output for up to 2 models makes a figure that reports
% spot price, futures prices and model-based futures price. 
 
textsize = 12;

yearticks = 2000:5:2020;
outputArg1 = figure(gcf);
colorimage = 0;
if colorimage == 1;
    firstcolor = [7,60,101]/255;
    secondcolor = [241,65,36]/255;
    thirdcolor = [94,204,243]/255;
    fourthcolor = [255,161,91]/255;
else
    firstcolor = [60,60,60]/255;
    secondcolor = [120,120,120]/255;
    thirdcolor = [150,150,150]/255;
    fourthcolor = [160,160,160]/255;
end;


    
    
cicolor = [0.85 0.85 0.85];


plot(xdates,spotprices,'-','linewidth',2,'color',firstcolor)  % Spot Price
hold on
plot(fdates,futures,'--','linewidth',2,'color',secondcolor)   % 2-yr ahead future
plot(xdates+30, modelfut,':','color',thirdcolor,'linewidth',3)  % 2-yr model implied
if nargin > 6
    plot(xdates+30,othermodelfut,'-.','color',fourthcolor,'linewidth',3)        % Particle Filter
end

if nargin > 5;
    lwbnd = [xdates+30,modelfutci(:,1)];
    upbnd = [xdates+30,modelfutci(:,2)];
    cipatch = [lwbnd;upbnd(end:-1:1,:)];
    ci=patch(cipatch(:,1),cipatch(:,2),[0.85 0.85 0.85]);
end

plot(xdates,spotprices,'-','linewidth',2,'color',firstcolor)               % Spot Price
plot(fdates,futures,'--','linewidth',2,'color',secondcolor)           % 2-yr ahead future
plot(xdates+30, modelfut,':','color',thirdcolor,'linewidth',3)  % 2-yr model implied
if nargin > 6
    plot(xdates+30,othermodelfut,'-.','color',fourthcolor,'linewidth',3)        % Particle Filter
end



vaxis = axis;
axis([xdates(1) xdates(end)+365 vaxis(3) 180])  

datetick('x','keeplimits');
if nargin < 7
leghandle = legend('Spot Price','2-year Ahead Futures Prices (Data)','2-year Ahead Futures Prices (Model)', 'Location','northeast')
else
leghandle = legend('Spot Price','2-year Ahead Futures Prices (Data)','Benchmark Model',legendtext, 'Location','northeast')
%legendtext choices could include 'Particle Filter Model' and Constant Gain
%Model \chi = 0.02
end

%orient landscape
ylabel('Dollars Per Barrel','Fontsize',textsize)
set(gca,'Fontsize',textsize)
set(leghandle,'Fontsize',textsize-3);
set(gca,'xtick',datenum(yearticks,6,30))
datetick('x','keeplimits','keepticks');
%orient landscape
%print -dpdf -bestfit newfigures/KF_wPF_2018



end

